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SUMMARY 


A recent application of numerical optimization to the design of advanced 
airfoils for transonic aircraft has shown that low-drag sections can be 
developed for a given design Mach number without an accompanying drag Increase 
at lower Mach numbers. This Is achieved by Imposing a constraint on the drag 
coefficient at an off-design M?ich number while minimizing the drag coefficient 
at the design Mach number. This multiple design-point numerical optimization 
has been Implemented with the use of airfoil shape functions which permit a 
wide range of attainable profiles during the optimization process. Analytical 
data for the starting airfoil shape, a single design-point optimized shape, 
and a double design-point optimized shape are presented. Experimental data 
obtained Iri the NASA Ames Two-by Two-Foot Wind Tunnel are also presented and 
discussed. 


INTRODUCTION 


The design of supercritical airfoils for advanced high speed aircraft has 
been facilitated by the computerized analytical methods which have been devel- 
oped In recent years. Although these methods provide good performance predic- 
tions for each individual design point which Is considered, they do not allow 
the designer to automatically consider off -design characteristics during the 
design process. A method which does provide multiple deiign-point capability 
Is described In reference 1 . It Is based on design by numerical optimization. 
An application of that method to a single design-point and a double design- 
point airfoil optimization Is addressed In the present study. The double 
design-point optimization produced a low drag supercritical airfoil 'or a given 
Mach number subject to a drag constraint at a lower Mach number. 

The treatment of the supercritical airfoil design problem by this method 
has been facilitated by the development of a set of airfoil shape functions 
(reference 1) which provide a wide range of attainable profiles during the 
design process. The coefficients of these shape functions are used as design 
variables In the numerical optimization technique which consists of two exis- 
ting computer codes: (a) an optimization program based on the method of 
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feasible directions (reference 2 ) and (b) an aerodynamic analysis program 
based on an iterative solution of the full potential equation for transonic 
flow (reference 3). 


SYMBOLS 


a^ shape function coefficients 

c chord 

Cjj section drag coefficient 

section lift coefficient 

C section pitching moment coefficient 

m 

Cp pressure coefficient 

f^ airfoil shape functions 

M Mach number 

X airfoil abscissa 

Y airfoil ordinate 


DESIGN METHOD 


Only a brief description of numerical optimization will be given here. 

A complete discussion of the technique can be found in reference 4. 

A schematic flow chart of the numerical optimization design program used 
during this study Is shown in figure 1. A baseline airfoil is required to 
start each design problem. The airfoil shape is represented in the program 
by the following equation: 


'' " 'basic * S 

where Ybasic the set of ordinates of the baseline airfoil and fi are the shape 
functions. The shape functions are added linearly to the baseline profile 
by the optimization program to achieve the desired design improvement. The 
contribution of each function is determined by the value of the coefficient, 
associated with that function. These coefficients are therefore the 
design variables. Other inputs to the program include Mach number, angle of 
attack, and any constraints to be imposed on the design. 
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. T*’®J’yPothetical design problem represented by the flow chart Is draa 

Jlli^a? wStherX^^^^ conJtraineJ to s^' lpecf?ied 

enother Jtech number* M 2 . The optimization program begins by chanainc 

JStuSr?!! *'*® Initial value of zero to 0.001 ^ It 

aerodynamics program for evaluation of the drao coefficient at 
both Mach numbers Mi and M 2 after each change. The valul S O.SSl ll sSLhJt 

JaHaSe? ?T *^ep changrin the deslST 

variables to calculate the required partial derivatives. The partial derl- 

^tives of drag with respect to each design variable form the gradient of drag 
y^ri* The direction in which the design variables are changed to reduce the 

active. The optimization program then Incranents the 
design variables in this direction until the drag starts to incre«e12cluse 

is eSwiStwed^ ^l/either^2f Mach number M 2 

IS encountered. If either of these possibilities occurs, new gradients are 

calculated and a nm direction Is found that will decrease drag without vio- 

yr ’ *•'“« of om for Mach nunber H? ll 

fin h^ biS icMiJiSI**' “ «2- ™<iu1re<l opt1«ized air- 


AIRFOIL SHAPE FUNCTIONS 


bv airfoil design by numerical optimization Is facilitated 

iLJrj® ® geometric shape functions, each of which affects a different 

limited region of the profile. General classes of such functionrwhich havf 

^0 Optimize supercritical airfoils are described In 
reference 1. The shape functions that were used In the present study were 

general functions and were applied to the airfoil upper 
exponential decay function and the sine functions are 
figure 2. The exponential decay function, fi, provided varia- 
tions in curvature near the airfoil leading edge. In the^^ne functions the 

perturbatlons^of f5°*'fr*f/SSH‘*f"®^®* assigned so that the maximum 

percuroaiions of f 2 , fa, f 4 and f 5 were at 20, 40, 60, and 80 percent of the 

chord respectively. The width of the region affected by each sine function 
was controlled by the localization power, 3. Previous studies (reference 1) 
have found that these shape functions provide a broad range of smooth airfoil 
contour modifications during the optimization process. 


ANALYTICAL DESIGN RESULTS AND DISCUSSION 


*1... importance of considering off-design performance of an airfoil durina 
the design process will be Illustrated by cSnpiirlng the results if a llngiir ® 

' '““I!’’’* '*«t*9'’-l><><"t optlmlaatlon. The first 
recontounng the upper surface of an existing supercritical airfoil 

if ® single design Mach number. The second consists 

?t thi nicini^M/h® “^E®'" sa">e airfoil to reduce the wave drag 

at the design Mach number subject to a drag constraint at a lower Mach number. 
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The calculated wave drag (reference 3) for Mach numbers near drag diver- 
gence for the starting airfoil and the two optimized airfoils are presented 
In figure 3. All these data are for 0,40 Ci, the design lift coefficient of 
the starting airfoil. Mach number 0.78 was arbitrarily selected as the primary 
design point, I.e., the Mach nundrer at which the drag would be minimized. 
Results of the single design point optimization are Indicated as 412M1. The 
drag at Mach number 0.78 Is significantly less than that of the starting air- 
foil and as a result the drag rise occurs at a higher Mach number. However, 
the drag at lower Mach nimbers, 0.76 and 0.77. Is greater than that of the 
starting airfoil. This local region of drag-creep could limit the usefulness 
of the Improved drag rise characteristics of the optimized airfoil. 

In order to avoid the drag-creep problem, the airfoil was optimized 
a second time with an upper bound of .0005 Imposed on the drag coefficient 
at Mach number 0.77. Results of this double design-point optimization are 
Indicated In figure 3 as 412M2. The drag rise for this airfoil occurs at 
a slightly lower Mach number than It does for 412M1, but there Is no drag- 
creep over the range of Mach numbers for which the airfoils were analyzed. 
Therefore, airfoil 412M2 is the more desirable design. 

Chordwise pressure distributions for the starting airfoil and for airfoil 
412M2 at Mach number 0.77 are presented In figure 4. The reason for the lower 
wave drag of the optimized airfoil Is obvious. The starting airfoil has a 
well developed shock at approximately 40 percent of the chord, but airfoil 
412M2 does not. Instead, It exhibits a gradual recompression from approxi- 
mately 10 percent to 50 percent of the chord. The geometric modification 
which has produced the pressure distribution change Is shovm In figure 5. 

This modification Is primarily a reduction In surface curvature from 5 per- 
cent to 40 percent of the chord. 

The aerodynamics code that was used in the optimization program Is an 
Inviscid, potential flow analysis method. In order to account for first 
order viscous effects 1n the flow field solution, a boundary layer displace- 
ment thickness was added to the starting profile before the optimization 
process. The displacement thickness was calculated for the pressure distri- 
bution of the starting airfoil at a Mach number near Its design condition, 

0.78. remained unchanged throughout the optimization process, and each 
of the optimized airfoils Included this same passive displacement thickness. 
Therefore, the analytical characteristics of the airfoils did not reflect 
potential changes In boundary layer behavior due to changes In the chordwise 
pressure distributions. 

Another aerodynamic analysis code (reference 5) was used to evaluate 
the active boundary layer characteristics of the starting airfoil and opti- 
mized airfoil 412M2. That computer program Is also based on an Iterative 
solution of the full potential equation for transonic flow, and it includes 
a momentum- integral calculation of the turbulent boundary layer parameters. 
During the solution, the airfoil geometry Is regularly updated with the 
boundary layer displacement thickness. The results of the viscous analyses 
with that code for Mach numbers between 0.76 and 0.81 Indicated that the 
differences In boundary layer cha’ uteri sties would be small. The calculated 
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wave drag for the starting airfoil and airfoil 412M2 is presented in figure 6. 
The relative increase in the drag rise Mach number is in good agreement with 
the results of the inviscid code (figure 3). 


EXPERIMENTAL RESULTS 


Models of the starting airfoil and airfoil 412M2 were tested in the NASA 
Ames Two-by Two-Foot Wind Tunnel. Data were obtained at angles of attack 
from -40 to stall at Mach numbers from 0,20 to 0.81. The test Reynolds 
number varied with Mach number as presented in Table I. Preliminary data 
from the test are presented in figures 7 and 8. The incremental values of 
drag coefficient, Cj,, have been referenced to the minimum drag measured for 
either of the airfoiis at each lift coefficient. Thereby, extraneous com- 
ponents in the absolute drag level have been excluded from the comparison. 

Drag characteristics for the starting airfoil and airfoil 412M2 at lift 
coefficient 0.40 (figure 7) indicate a difference in drag ri'je Mach number 
of 0.02 to 0.03. This improvement is greater than had been predicted by the 
analytical codes (figures 3 and 6). Drag characteristics for the two airfoils 
at lift coefficient 0.60 (figure 8) also ^indicate significantly less drag for 
the optimized airfoil 412M2 at all Mach numbers* Therefore the airfoil per- 
fonnance at this off-design condition has not been adversely affected by the 
design Improvement at 0.40 lift coefficient. 

The low speed drag-creep which occurs between Mach numbers 0.60 and 0.70 
for both airfoils is caused by the initial development of supercritical veloc- 
ities over the upper surface and the formation of a mild shock n<»dr the leading 
edge. Only Mach numbers greater than 0.76 were considered during the present 
analytical design study, but the numerical optimization technique could also 
be applied to the minimization of drag-creep at the lower speeds. 


CONCLUDING REMARKS 


A technique for designing low-drag supercritical airfoils has been demon- 
strated. The technique was used to modify the upper surface of an existing 
12 percent thick supercritical section to achieve a substantial drag reduction 
at Mach number 0.78 without an accompanying drag Increase at lower Mach numbers. 
The ability to treat this and other multiple design-point problems has been 
achieved by the use of a set of airfoil shape functions which provide the 
necessary flexibility in the profiles that are attainable »/y the numerical 
optimization design technique. Such capability is Important because each 
design point might require the modification of a different region of the profile. 

The two design-point problem considered in the present study illustrates 
the advantage of design by numerical optimization. Aerodynamic requirements at 
any number of off -design conditions are handled automatically without manual 
Intervention by the designer. Therefore, it provides a powerful tool for the 
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design of airfoils to meet specified performance goals throughout a flight 
envelope. 
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TABLE I. 

Mach Number 

0.20 

0.40 


- MIND TUNNEL TEST PARAMETERS 
Angle of Attack 

BS!!S£ Reynolds Number 

-4° to stall 1.9x10® 

3.0 X 10® 

4.0 X 10® 


0.60 



aerodynamic program 


OPTIMIZATION PROGRAM 


INPUT INITIAL AIRFOIL AND 
DRAG CONSTRAINT AT M 2 

initialize design VARIABLES 

(ai - . . a„> 


CALCULATE Cj AT Mj AND Mj 
CALCULATE C^ AT Mj AND M 2 ■ 



, PERTURB DESIGN VARIABLES 




\ *Cd 

CALCULATE-^' 

FOR Ml AND M 2 

FORM GrIoIENTS VCjATMj 
AND 7CdATM2 
.SET Oj.*. 

\/AI IIES AS determined BY GRADIEIffS 


On TO NEW 


CqATMi IS MIN. 

Cq a/m2 IS LESS 
THAN SPECIFIED 
VALUE 


pji 


CqATMjIS 
NOT MIN. 

AND/OR — 
Cj)ATM2 is 
GREATER TWN,^ 
specified VALUE 


OUTPUT FINAL 
PROFILE AND 


Figure 1.- Flow chart of numerical optimization design program. 



Figure 2.- Airfoil shape functions. 
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Figure 5.- Airfoil geometry comparison. 


COMPUTER CODE: GARABEOIAN-KORN WITH BOUNDARY LAYER 



Figure 
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Comparison of wave drag 

airfoil and optimized airfoil A12MZ. 


for the starting 
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